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Abstract 

A statistical field theory for an inhomogeneous liquid, a planar liquid/vapor interface, is de- 
vised from first principles. The grand canonical partition function is represented via a Hubbard- 
Stratonovitch tranformation leading, close to the critical point, to the usual cj)'^ scalar field theory 
which is then rigorously considered at the one-loop level. When further simplified it yields the 
well-known capillary wave theory without any ad hoc phenomenological parameter. Internal coher- 
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numerical simulations is stressed. 
PACS: s61.30.Hn; 64.70.F; 68.03. Cd. 

Keywords: Field theory; Gas-liquid interface; Surface tension. 

* russier@icmpe.cnrs.fr 

^ jean-michel.caillolQth.u-psud.fr 



1 



1 Introduction 



The structure in the interface between two fluid phases at coexistence plays a key role in many specific 
situations, such as, for instance, the wetting transition and related phenomena [1-7]. An understanding 
and proper description of this structure is of course challenging and presents a great interest in itself 
from a theoretical point of view. For simple fluids, one of the key features of interfaces that separate 
subcritical phases at coexistence is the presence of thermally activated capillary waves [8] coupled 
with density fluctuations in the bulk. Since the pioneering work of Buff et al. [8] capillary waves have 
been widely studied and are usually described through the introduction of an ad hoc effective surface 
Hamiltonian, written as a functional of the interface height, Zint{s), s being the coordinate parallel to 
the interface -in the case of a planar interface as considered henceforth-. H[zint{s)\ can be obtained 
through different ways and is most generally deduced from phenomenological arguments taking into 
account the interplay between gravity and surface tension effects. It is noteworthy that most works in 
this area are devoted to the wetting transitions where the Hamiltonian includes an external potential, 
say W{zint{s))-, which pins the interface to a solid substrate. 

The Hamiltonian H[zint{s)\, can be treated in the framework of statistical field theoretical meth- 
ods [2,4], including renormalization group (RG) [9] technics. For instance, recently, the hierarchical 
reference theory (HRT) has been generalized to inhomogeneous cases in order to deal with the wetting 
transition [7]. In the case of the liquid/vapor interface, when H[zint 

(s)] is treated at the gaussian 

level, the capillary wave theory (CWT) is recovered, leading to the well known behavior for the 
height-height structure factor SazAz{q)- Moreover, 3D cj)'^ model in which the stochastic variable <j) 
does not coincide with the height Zint{s) were studied in refs [5,6,10,11]. Alternative ways to deduce 
an effective Hamiltonian from a density functional theory (DFT) including the effect of gravity [12] or 
the local curvature of the interface [13, 14] have been also considered. Recently [15], a model for the 
density profile based on an extension of a displaced profile approximation, i.e. in which the profile 
is written as a function of [z — Zint{s)), was considered. In addition, the model takes into account, 
together with the surface fluctuations described by the height Zint{s), bulk phase fluctuations. 

The purpose of the present work is to provide a simple description of the liquid / vapor interface 
of fluids in the framework of an exact statistical field theory. The latter is obtained from a Kac- 
Siegert-Stratonovich-Hubbard-Edwards (KSSHE) [16-20] transformation devised to rewrite the grand 
canonical partition function (GCPF) of the inhomogeneous fluid as the partition function of a statistical 
field theory involving the stochastic real field cp. The method was mainly used for homogeneous fluids 
in refs [22,23] and is considered here to study a planar interface. The theory is studied at the one-loop 
level where everything can be achieved analytically. The link between correlations of the field and 
density correlations are established as well as the expressions for the surface tension and the surface 
structure factor. Technically, this is done via the determination of the eigenvalues and eigenfunctions of 
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the second functional derivative of the mean field KSSHE Hamiltonian, which can be done exactly, as 
well known, for the one- dimensional kink [21]. Taking into account the whole spectrum of eigenstates 
proves to be of the outmost importance to obtain the correct result. We thus obtain a satisfying 
picture of the liquid vapor interface, at least qualitatively, even at this simple level of description. One 
of the salient points of our work is that we recovers both the CWT and its first extension, namely 
the appearance of the bending rigidity factor k and the coupling between surface and bulk density 
correlations in S{q), without invoking any ad hoc phenomenology. 

2 KSSHE transform and mean- field approximation 

The important steps leading to the appropriate statistical field theory for a liquid are outlined in 
Appendix [Aj the interested reader is referred to ref. [17,22,23] for deeper details and to ref. [24] for 
an alternative formulation. We consider a simple fluid whose pair interaction potential includes a hard 
sphere (HS) repulsive part and a soft attractive part, denoted by v(r). Let w{r) = — /3v(r) as usual 
and suppose that only purely attractive potentials are considered (w is a positive definite operator, 
i.e. w{k) > 0). We work in the grand canonical ensemble, namely at constant chemical potential 
u = and we consider the grand canonical partition function (GCPF) related to the grand potential 
S = exp(— Following the notations of appendix lAl we obtain for H 



where (j) is & real scalar field, Vcp a functional integration measure, v = v — w(0)/2 and, finally, Hhs[i^] 
is the hard sphere grand partition functional of the local chemical potential v(f). A Landau-Ginzburg 
form is obtained from a functional Taylor expansion of In(HHs) around a conveniently chosen reference 
chemical potential of the hard sphere fluid, (see Appendix [Aj. The propagator is expanded up to 




(1) 



order k? 



in Fourier space and we are left with 



I 




(2a) 




(2b) 



n>3 




(2c) 
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where the couphng constants Kq,K2 and u„ depend only on the equation of state of the hard sphere 
fluid and the potential w (see below). This formulation allows an exact mapping between the densities 
and their correlations on the one hand and the mean value and the correlations of the field on the 
other hand, as examplified in Eq. ()A-13p in the appendix. For instance at coexistence, the densities pi 
and pg of the two phases correspond to the two values < 4> >i and < 4> >g. The value of uq is chosen 
in such a way that u^{i'o) vanishes and the coexistence condition between the liquid and vapor phases 
in the absence of external field is = 0. Kq is related to the deviation from the critical temperature: 
t = {T,-T)/T, cc i-Ko). 

In the inhomogeneous system, we assume a special realization of the two phases coexistence: we impose 
explicitly the occurrence of a bulk liquid and a bulk gas, at densities pi and pg separated by a planar 
surface located at z = zq. The system is supposed to be bounded in the direction z according to 
{zq — L) < z < (zq + L), L being macroscopically large though formally finite. The localization of 
the interface is obtained through the introduction of a small valued external field hext{z — zq), odd in 
{z — Zq). Since the mean field profile takes constant values in both bulk phases, so is the case for hext 
which therefore plays the same role as the truncated gravitational field used in [12]. The advantage 
of our choice will appear below; we emphasize that we do not aim to study gravitational effects but 
rather to fix the location of the interface at z = zq. The inhomogeneous mean field theory or saddle 
point equation, 6H[(p]/64> = hext{z — zq), is solved by exploiting the fact that we look for a monotonous 
solution for the mean field profile 4>c{z) = 4'MF{z)/(f)b, where (pb is the value of (/>mf(-z) when z zq + L, 
and by expanding h^xt as a function of (l)c{z), namely hext{z) = hext{<Pc{z)) . The mean field equation 
then reads 



K2^-^ = Ko(f>UF + ^{(t>MF) - hext{(pc) (3) 



Since we consider a local approximation for the n > 3 HS kernels (see Eq. (|A-14p ). the r.h.s of Eq.([3|) is 
a polynomial in (pMpiz — zq). Limiting V{(f)) to the first non vanishing term, namely V{cj)) = (n4/4!)(/)^, 
the analytical form for (/)mf is not modified by the introduction of an external field if the later expands 
as h{(pc) = hi4>c + (/i3/3!)</>c (/ii, /13 constants), as already noted by Zittartz [25] when h{4>c) = /Ji</>c 
(This results holds in a more general way for a potential V{4)) even in 4> and an external field odd in 
(f)c where the highest power of V{(j)) is that of /i((/'c) plus one). The solution of the mean field equation 
in presence of the external field, (/>MF(-f^2> -f^Oi hext), then coincides with the solution in the absence 
of external field, with modified values of the coupling constants Kq and M4. 

with 6K0 = /ii and 6u, = hs (4) 
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The solution of equation ^ is 4'mf{z) = (j)b4'ci{z — zq)/1) with (pdx) = tanh(x), = {—6Kq/u'^)^^'^ 
and / = {2K2/ — K'qY/"^ is the bulk correlation length which coincides also with the intrinsic interface 
width. We emphasize that we deal only with subcritical temperatures where K'q < 0. We first have to 
check the mapping between the physical external field seen by the liquid, say ^(z) = PVext{z), and the 
external field which enters the effective Hamiltonian of equation (f2al) . We thus have to solve 



(5) 



B(l) = B„-/„r.(i,2)*„.(2M2.,.,,(l) 



Bq given by Eg. ([2c]) is the constant external field appearing in the effective Hamiltonian when hext 
vanishes. The coexistence condition at constant chemical potential still reads Bq = 0. Since the 
propagator has been expanded up to order k'^ (see appendix [A|, is given in real space, at the same 
order of approximation, by 



1 

Wq 



1 + 



W2 
Wo' 



(6) 



where wq and 1112 denote the coefficients of the expansion at order k'^ of the Fourier transform, w{k), 
of w. By limiting ourselves to functions depending only upon z, we thus have 



J Wo 



2 1 



1 + 



W2 d 

Wo dzf 



"^extil) = hext{zi) 



(7) 



We introduce the bulk correlation length, I and we look for a solution depending on t = (pdz/l) for 
both hext{zi) and ^ext{zi)- Equation ([7]) then reads 



2\2 



5t2 



d , ,2^0 



dt W2 



W2 



I h^xt 



(8) 



A one to one mapping with hext limited to the third power in (pc is obtained when d^'^/dt^ = and 
we get 



-^{z) = ^'l(/>c(z) ; hextiz) 

hi 
h. 



Wo V ^0 

2W2 
P Wq 



(9) 



Formally I is finite but clearly the present approach should be valid only but in the vicinity of the 
critical point where / 00 and we can therefore neglect the term /13. 

The mean field surface tension follows from the identity ilMpi'^] = H[(j)MFiz)] which is easily calculated 
by using both the equation satisfied by cpMF and the expression of the external field contribution which 
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takes a simple form with the particular choice for hext, with the result 

1 /2 

Pa^MF = Aa\-2K,K'i/uifl^ + 2a' h, {l + 

which is similar to the expression given by Brilliantov [11] in the absence of the external field and 
differs from that given by Zittartz [25] because of a different choice for the external field as a function 
of (j)c. When /i3 = and /ii — > at finite Kq we get 

l3a'-iMF{Ko, K2, hi^O)^ Pa'jMFiKo, K2, h = 0) + (10) 

For i — > 0, one recovers the mean field exponent, ^mf ~ t^^'^- 



3 One loop equations 

In order to go beyond the mean-field approximation, we expand the Hamiltonian about i?i>MF(^) and 
consider the Gaussian approximation. Prom here, we drop the cubic term in the external field {i.e. we 
set /13 = 0) and, in order to unclutter notations, we define hi = hi/cph. Without any lost in generality, 
we also choose zq = 0. The first correction stems from the second order term and we have [zq = 0), 

H[cj) = ^UF + x] ^ H[<Pmf] + IJ x{l)H^^Hl,2)x{2)dld2 (11) 

The operator H^'^^ = 5'^ H / 5(j){l)5(j){2) is diagonalized, after a Fourier transform parallel to the surface, 
in the set of the eigenfunctions ipx solution of the eigenvalues equation (e;^ = K2{q^ + hi + uJx/l'^)) 

i^2(-^ + 6(tanh2(z*)-l) + (4-a;A))^A(^*) = z* = z/l (12) 

the solutions of which are known [5,25,26]. We emphasize that keeping only the linear term in the 
external field leads to the same Schrodinger - like equation ([T2]) as in the absence of external field 
and the only change is a shift of the eigenvalues e\ of hi. The spectrum of eigenfunctions includes 
two bound states, ifQ = Co/ cosh^(z//) and ifi = Ci sinh(z//)/cosh^(z//), with ojq = 0, u>i = 3 
respectively and a subset of unbounded states, or continuum spectrum ^k{z), with kl = ^lj^. — 4 which 
behave as plane waves in the bulk phases, i.e. far from the interface, and are given by 

ifkiz) = C{k)e'^^ [2 - k'f - 3i kl tanh(z//) + 3 (tanh2(z/0 - 1)] (13) 
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Notice that the shift of the eigenvalues due to the external field is of crucial importance especially for 
the lowest bounded state (^o the energy of which eo tends to a finite value, ^2^1 when g — > 0. As we 
shall see in the sequel this corresponds to a pinning of the interface by the field. It is important to 
note that ^po{z) oc d(j)MFiz)/dz. The constants Co, Ci and C{k) are determined in order to normalize 
the ipx on the interval [—L,L]. The whole spectrum of eigenstates must then be orthogonalized. For 
the bound states this is a direct consequence of the eigenvalues equation (fT2]) in the limit L/l — > 00, 
i.e. disregarding terms or order exp(— 2L//). On the other hand the orthogonalization of the subset 
of unbounded eigenstates is no more a consequence Eq. (fT^ but rather follows from the boundary 
condition 

[VkVl'tL = for k^k' (14) 

from which we get the non trivial dispersion relation leading to the density of states [x = k I) 

L 



n{k) = 

vr TT 



The closure relation then follows 



+ 



(l+x2) (4 + x 



+ lf{x) (15) 
vr 



^ ^ [n{k) \^k{z)\^ - Ij + ifliz) + ^l{z) = (16) 

It is important to note that f{x)dx = -2 which proved useful in all calculations. Zittartz [25] have 
already got this density of states but without explicitly mentioning the need of the orthogonalization. 
In [27] a similar kind of dispersion relation was obtained in the modeling of the charge density profile 
of electrolytes in the framework of another field theory. We emphasize that the orthogonalization of 
the eigenstates is of crucial importance for the calculation of gaussian functional integrals. In order to 
calculate ln(H) and the correlation functions, we write the fields xi^ in the basis which diagonalizes 
the operator H^'^\ namely a Fourier transform in the s- direction and a projection on the (px{z) 
where {A} = {n = 0, 1; k} denotes the whole spectrum of eigenstates. From usual gaussian functional 
integrals [28] we get 

ln(H) = -H[cPmf] + ^ / In — (17) 



2 J (2^)3 V ^kiq) 
-f /(||f [E (<?))+ / dkfikl) ln(6,(g)) 



where the coincide with the eigenvalues of the second functional derivative of the Hamiltonian for 
the (j)^ homogeneous model taken at (j) = (pupi^ L). Therefore, the volume term of ln(H) in (fT7|) 
is nothing but (3pV at the one loop approximation [22]. We finally write the surface term of ln(H) in 
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term of the surface tension 7 = [ln(Hf,) — ln(E)]/S = ^mf + 7''^^- We obtain (7^^*-' = /37^^V^) 



v{l*) 



1 

8^ 



{t - l)ln(l 



^^t + l 



+ (t - 4) ln(l 



where qm and are the lower and upper bound respectively of the integral over q and a = (4 + /ii). 
A similar result was obtained in [25] for a spin model. This result differs from that obtained in [4] 
where only the n = eigenstate is kept. In the limit hi 0, we get 



7(i*)(/ii) ^ -/^^*\hi = 0) - {l/8TT)hi ln{hi). 



(19) 



We now consider the calculation of the field two-body correlation functions, 
G(^(l,2) =< x(l)x(2) >H- More precisely, we focus on the Fourier transform parallel to the surface, 
G(j,{zi, Z2,q)- For this we have to calculate the sum 



nt^i 'n{q) J e,{q) 



A contribution of the integral over k cancels exactly the direct contribution of the two bound states, 
which shows, once again, that the approximation consisting in keeping only the ipo{z) eigenmode is not 
sufficient. The result is (2:>,< = sup,inf{zi, Z2)) 



G^{zi,Z2, I 



91 exp(- |zi2| (x2 + 4)1/2/0 



2K2 2;2(x2 + 3)(2;2 + 4)1/2 
1 + x2/3 + (x2 + 4)1/2 tanh(z>//) + tanh2(^>//) 
1 + a;2/3 - (x2 + 4)1/2 tanh(z<//) + tanh2(z<//) 



(20) 



where x'^ = {q^ + hi)P. Notice that formally this expression coincides with that obtained in the absence 
of the external field, the only difference being the definition of x which takes a finite value, hi when 



4 Results and discussion 

4.1 Capillary behavior and surface structure factor 

We started from the expansion of the effective Hamiltonian on the basis of the eigenstates ^\{z). The 
first eigenstate, (po{z) is proportional to the derivative of the mean field result (pc{z), the proportionality 
constant being determined by the normalization. If we keep only this first eigenstate, the expansion of 
X reads x{s, z) = £^{s)ipQ{z). Hence the field takes the form 4){z) = 4>c{z — Zint{s)) and only the linear 
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term in the expansion of (j) with respect to Zint ■ remains. This corresponds indeed to the so-called 
rigidly displaced profile approximation, where Zint = —o,^,{s) represents the fluctuating location of the 
interface. The corresponding contribution to H is a functional of Zint 

{s) which defines an effective 

surface Hamiltonian given by 



where we have used eo = K2{q^ + hi) since = 0. This Hamiltonian pertains to the capillary wave 
models [12,29-31]. In general the capillary wave Hamiltonian is obtained either from phenomenological 
arguments [29] or by restricting the fluctuations beyond the mean field level to those deduced from a 
rigidly displaced mean field profile [30]. In [12] a formulation based upon the direct correlation function 
of the liquid state theory is used. In [32] a supplementary Zint{sY term in the surface Hamiltonian 
([^Tl) was introduced in order to account for curvature effects of the interface; this additional term must 
vanish with the external field, [31] and in any case cannot be generated in the framework a one loop 
scheme. From the normalization of tpQ and the analytic form of the mean field profile we deduce K2/a'^ 
=(t'^Pjmf — o''^l4>b^i and accordingly we rewrite ([2T]) in the form 



which coincides exactly with the usual effective surface Hamiltonian iJcwT[-2mt(s)] of the CWT the- 
ory for the free surface, when the external field vanishes. Therefore, in the limit /ii — > we obtain 
the CWT as the lowest approximation beyond the mean field approximation without invoking phe- 
nomenological arguments. The structure is characterized by the height-height correlation function, 
< Zint{si) Zint{s2) > 01 its Fourier transform parallel to the surface which defines the surface struc- 
ture factor, SAzAzig), where Zint{s) is the location of the interface relative to its mean value. We 
consider / xis,z)dz/A(l)f,, where Acpb = 0mf(-^) — </'mf(— -^), as a measure of the instantaneous lo- 
cation of the surface at s, which amounts to define the location of the surface from a constraint on 
the integral of Xi as is done in [15] for the density profile. We are then led to identify S^zAzil) 
= (Acpi,)"'^ J G(p{zi, Z2,q)dzidz2, which corresponds to the Sic used in [15]. It is important to no- 
tice that the coupling with the bulk fluctuations are included in the present formulation through the 
eigenstates of the continuum. We also introduce the effective surface width, or surface corrugation, 
a'^-^f = < Zint{si) Zint{si) >. Hcnccforth we consider only small values of hi and more precisely we 
relate the value of hi to an effective lateral size of system, according to hi ~ as it will be justified 
from the analysis of the capillary wave limit. The behavior of Sa.z Az{q) is analyzed from the function 




(21) 



= I [PlMF - a^l^hl] I ds[{ds{Zint{m^ + hlZ^nt{^^] 



(22) 
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g{q) = J dzi J dz2G^{zi, Z2, q). Prom ([20]) we get for the leading term of G^{zi, Z2^q) when g — > 

K2[q^ + hi) lK2q^ + hi 

(23) 

The first term of the r.h.s of Eq ([23]) corresponds to the CWT hmit especially when hi = 0, and the 
correction, which vanishes at hi = 0, plays a role only far from the interface where ipQ{zi)ipo{z2) is 
negligible and moreover when the two particles are on the same side of the interface, i.e. when e"^'^^^'/^ 
is not negligible. Thus we can rewrite the g — > behavior of G^{zi, Z2,q) as 

G^{zi,Z2,q^0) ^ ;^^(^^^'^o(^i)'^o(^2) + ^/e-2|--l/^ (24) 

where / takes a constant value, / = 1, when both jzil , |2;2| >> Z and are of the same sign. The second 
term of the r.h.s of ([24l) is nothing but the limit at g — > of the bulk correlation function and is not to 
be included in the capillary limit of the model. In any case, for sufficiently small values of hi, the bulk 
contribution to Eq. (p^ becomes negligible at g = 0. This is made more precise below concerning the 
behavior of SazAz- Let us denote by Gcw the first term of the r.h.s. of Eq. (f24|) and by gcwiQ) its 
integral over zi and Z2. It is easy to show from the relation already used between the normalization of 
(po{z) and 7a/f that 



{AcP,y'^'' ^ (7MF - crH^t^lhiW + hi) 

In the limit /ii = this is exactly the CWT behavior, leading to the well known logarithmic divergence 
of the squared surface corrugation for which one gets (cr^-^-^)^ = {At^'^mf)~^^^{<i\j/q'^)- The effect of 
the external field is to make the integral over q entering in the determination of a'^^f finite at g ^ 
with the result 

K^^)^ = ^2 1- [1 + 0{a'hi)] (26) 

where we have assimilated qM to the inverse of the molecular size, 2iT/a. Since the result for the CW 
at the free interface in the absence of external field is a logarithmic divergence with the lateral size, 
~ ln(cr/La;) [8], this shows that, as we already mentioned, the effect of the external field, is to pin 
the interface in such a way that the capillary wave are restricted to a lateral scale ~ ''^^^^ . It 
can be shown from pOl) that the contribution to g{q) diverging with the system size, due to the bulk 
correlations, is exactly (2L) times the integral of 0^(^12, g),the correlation function of the bulk phase. 
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over zi2 If we keep only these two terms we get 



Saz Az[q) ^ -r-r2 / Gb{q, Zujdz^ + ( 2 , T. \ + '^^^l) 

A* J 7A/F(r + /ii) 



Given that the bulk term leads to a constant when g — > 0, we see that (l27l) presents a cross-over 
like behavior in terms of wave vector q, where a threshold value Qs naturally appears, separating the 
capillary wave behavior at small values of q from the bulk like behavior at q > Qs, with qg given by 

We see that this cross-over behavior holds only if hi is negligible compared to the r.h.s.of Eq. ([28]l or 

in other words 

< 2L P^mfP 

and confirms the above mentioned separation of the capillary behavior from the bulk correlations when 

~l/2 

this condition is fulfilled. We emphasize that can be understood as the inverse of an effective 



lateral size pinning the interface. The cross-over like behavior resulting from (|28p is in agreement with 
that of Refs. [15,33,34]; however, in the present formulation the bulk fluctuations come out of the 
calculation through the continuum subset of eigenstates and have not to be added to the field profile. 
Furthermore, from (|27l) we can drop exactly the bulk contribution, and doing this we define a purely 
interfacial contribution to Sa.zAz{(i) (see fig.([TJ). 

Sintiq) = Saz Az{q) - J Gk{q, Zi2)dzi2 (29) 

Then the departure of Sint{q) from its 1/g^ behavior allows us to isolate the deviations from the 
capillary wave like behavior of SAzAz{q)- At this point, since we focus on the deviation form the 
capillary - like behavior, we can consider hi = 0. We define I3{q) = {q^ Sint{q))~^ —JMF- I3{q)/q^ is 
found nearly constant and thus lead us to define a bending rigidity of the interface, 

^ = x.M{i/q^)m) ^ {{Vq^)m) (30) 

which takes a positive value, as it should be for the stability of the interface [14], as is also found from 
the simulation results of Ref. [34] after a separation between bulk density correlations and interface 
fluctuations. In other words, the small q behavior of Sint{q) is 

Sint{q) ^ l/[7AfFg^ + K-q^] 
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Here we find n = k*(7mf<7^)/((7c)^, with the reduced value k* = 0.288. 



4.2 Surface tension 

In this section we focus on the behavior of the one loop contribution to the surface tension 7 in the 
case hi = 0, since its dependence on hi then vanishes (see Eq. (fT9]) ). The contribution 7^^^, given by 
p8]) . depends on the two bounds qm and qm of the wave vector q which are related to the relevant 
parameters of the interface: on the one hand IqM = ^nl/a where l/a \s the intrinsic width of the 
interface in unit a, IqM G [l)C>o[. On the other hand, Iqm = (2vr//La; where is the system actual 
lateral size; hence Iqm = {lqM){<y / Lx)- Thus IqM appears as a natural parameter, with IqM G [27r, 00] 
and l/{lqM) oc \/t. We can now re-write 7^^*) is a more convenient form {a / Lx G ]0, 1]) 



2{lqMf 



l{^+{lqMf)-l{'^+{lqM-r) 



(31) 



(Lx/a) is either the actual lateral system size or the scale at which 7 is measured, for instance in a 
numerical simulation (see ref. [35]), but in any case does not depend on t. We can note that whatever 
the values of a/Lx or (IqM), 17^7^^^ remains finite and more precisely 

■ri':ljkM - ~) ~ (32) 

which when t — > as t^^^. The results for a'^j^^*^ are displayed in fig. We interpret 

a'^^i^*) (^a/Lx) as the g— dependent contribution to 7 with q = (2-k(t/Lx). This means that the flat 
interface corresponding to the mean field approximation is obtained when no fluctuation at all are 
taken into account, namely for q ~ qM- This differs from what is done in Ref. [14] (see also ref. [36]) 
where the contribution to 7 due to the surface fluctuations vanishes at g — > 0. The small q behavior of 
7^^*) is easily obtained from (l3T]l and yields 



7S-7^^*^(0) + 77-[(9f^)^[a + ln((T//)]-2(g(j)2ln(gCT)] with a = 3.48491 



(33) 



from which we can deduce a crossover value of q given by qocr = {a/l)e°'/'^ ^ separating the q^ behavior 
from the (go")^ ln(g(T) behavior obtained for q < qo and q > qo respectively. It is important to note 
first that we always get an increasing 'y^^\q) and secondly that since qo is proportional to we get 
a plateau (corresponding to the q'^ dependence) at small values of a/Lx only when IqM takes small 
values, i.e. for the lowest temperatures (see figure ([2]). This behavior is in qualitative agreement 
with the simulation results of [35]. The term proportional to q^ in the variation of 7*^^-' with q can be 
interpreted as resulting from the energy necessary to bend the interface, and should be related to k 
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obtained from the behavior of Sint{q) (see eq. ([30]) . This is not a priori the case since we do not expect 
a fully coherence in the framework of a loop expansion between the energetic and structure quantities. 



4.3 Density profile and density correlation function 

Now we go back to the mapping between the averaged field, (j), and the density in the fluid, p on the 
one hand and between the field and the two body density correlation function G^"^ on the other hand. 
Formally we have from Eq. ljA-lSP ) 

pil) = pnsWo]il) - hextil) + J w-\l,2)<l)MFi2)d2 
0^^X1,2) = j j w-\l,l')w-\2,2')G^{l\2')dl'd2' -w-\l,2) (34) 

We recall that Gf'^\l,2) is p{l)p{2) h{l,2) + p{l)5{l,2) where h{l,2) is the usual total correlation 
function of liquid state theory [37]. We easily get 

p(l) = /OHs(i^o) + ^PiMzi) + ^pMzif (35) 



where 



Ap3 



1 



W2 



(36) 



If we limit the expansion to the term linear in hi we recover the usual tanh - like profile {pi + Pg)/2 + 
tanh(z — zq){pi — Pg)/2. Concerning the two body correlation function, since we limit the effective 
Hamiltonian to a (f)'^ model, the bulk field correlation function is the same in the two phases and as 
a result of equation ([3l|) this is also the case for the density correlation function Gp (1,2). This 
is of course a drawback of the simplicity of the present model for the description of the liquid vapor 
interface. To overcome this drawback, we may have to consider a true inhomogeneous reference fluid, 
which impose to chose two different values for the hard sphere reference chemical potential, uq. This 
is beyond the scope of the present work. 



5 Conclusion 

We have shown in this article that the one-loop approximation of the KSSHE statistical field theory 
of simple liquids gives a coherent and qualitatively exact description of a planar liquid/vapor inter- 
face. The mathematical treatment of the kink is subtle and requires to take into account the whole 
spectrum of eigenvalues and eigenstates of the Gaussian Hamiltonian and not, as often and incorrectly 
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proposed in the literature, to restrict oneself to its fundamental state. Moreover, without invoking any 
phenomenological description of the interface, we were able to recover the usual CWT surface Hamil- 
tonian as an approximation of the one-loop calculation. We also show that we can extract a rigidity 
bending factor which takes a positive value, in agreement with the requirement for the stability of the 
interface with respect to fluctuations. On the other hand, we clearly control the limit of the model by 
computing the two-body density correlation functions which appear to coincide in the liquid and vapor 
bulk phases, a drawback due to the identity of the two phases in the field (j) representation. 

Nevertheless, qualitatively, our results are in good agreement with numerical simulations. We 
stress that no parameter is involved in the theory which can be worked out for any pair potential. A 
quantitative agreement is expected only for long range pair potentials which we plan to investigate by 
numerical simulations in future work. 



A KSSHE Transform 

In this appendix we review the so-called Kac-Siegert-Stratonovich-Hubbard-Edwards (KSSHE) [16-20] 
transformation devised to rewrite the GCPF of a classical fluid as the partition function of a bosonic 
statistical field theory; more details are given in refs [22,23]. 

We consider the case of a simple three dimensional fluid made of identical hard spheres of diameter 
a with additional isotropic pair interactions v{rij) {rij = \\xi — Xj\\, Xi position of particle "i"). Since 
v{r) is arbitrary in the core, i.e. for r < fi, we assume that v{r) has been regularized in such a way that 
its Fourier transform v{q) is a well behaved function of q and that v{0) is a finite quantity. We denote 
by C M'^ the domain of volume V occupied by the molecules of the fluid. The fluid is at equilibrium 
in the GC ensemble and we denote by /3 = 1/kT the inverse temperature [k Boltzmann constant) and 
H the chemical potential. For the sake of generality the particles are subject to an external potential 
tp{x) and z/(x) = f3{fi — ^(x)) is the dimensionless local chemical potential. We stick to notations 
usually adopted in standard textbooks on the theory of liquids (see e.g. ref. ( [37])) and denote by 
■w{r) = —j3v{r) the negative of the dimensionless pair interaction. Moreover we restrict ourselves to 
attractive interactions, i.e. such that their Fourier transform w{q) > is definite positive for all q. 

In a given configuration C = (A^; xi . . .xn) of the GC ensemble the microscopic density of particles 
reads 

N 

p{x\e) = Y,5^^\x-Xi) , (A-1) 
and the GC partition function H [v] can thus be written as 

"M = £4t/ dl...dn e^Y>{-PVns{C) + ]-p-wp + V-i] , (A-2) 
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where p ■ w ■ p is a short hand notation for / p(l)tt;(l, 2)/9(2)(il(i2, i = Xi and di = d^Xi. For a 
given volume V, H [u] is a function of f3 and a convex functional of the local chemical potential u(x) 
[38,39]. In eq. (IA-2P exp{— (3Vi{s{C)) denotes the hard sphere contribution to the Boltzmann factor in 
a configuration C and V = u + vs where us = —w{0)/2 is (3 times the self-energy of particles. From 
our hypothesis on 'w{r), us is a finite quantity which depends however on the regularization of the 
potential in the core. 

We now recognize the Gaussian integral [40,41] 

exp • w-p^=^ jvfexp (^-^f ■ ■ ip + p ■ ip^ , (A-3) 

where ip is a real scalar field, Vp) a functional measure, w~^(l, 2) is the inverse of w{l, 2) in the sense 
of operators and the the normalization Mw reads as 



j Vp exp {^i^'f -w ^ -p^ 



V r d^k 
2 J (2^ 



exp ( -- / j:^lnwik) ] , (A-4) 



It follows then that the GCPF H [u] can be re-expressed as the functional integral 

E[u] = jvp expi-H[p,u]) (A-5) 

where the KSSHE action (or effective Hamiltonian) is expressed as 

H[p,u] = ^p-w~'^ -p -In {Ens [^ + ^]} ■ (A-6) 

The action H [p, v\ is non-canonical in the sense that the coupling between the field p and the external 
source v is non-linear as in usual Landau- Ginzburg actions [41]. In order to get a canonical field theory 
we perform the translation = -|- Az^ with Az/ = v — vq where i^o is some arbitrary uniform reference 
chemical potential to be chosen conveniently later. The Jacobian of the transformation is obviously 
equal to one and then, performing a Taylor functional expansion of the grand potential In {Hjjs [^^o + '^]} 
about vi^ we obtain a standard scalar field theory [22,23], characterized by the partition function 

E^ {B\ = AC' j V(p exp 00 ■ A-i - y [</)]+ i? • 0^ , (A-7) 
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where the propagator A{x,y) and the interaction V[(j)] are given by 



A-\x,y)=w-\x,y)+G'il^ii^o\x,y) , (A-8a) 

^M=-X]— dxi . . .dXn G^^l^{uo\xi, . . . ,Xn)cl){xi) . . .(j){Xn) , (A-8b) 
n=3 

where the ^ are the truncated n-body correlation functions of a HS fluid at chemical potential 
1^0 [37] 

.(n.)T^^, i^^ ^ ^ _ In Shs Wo 

6uoixi) . . . duoiXn) ' 



Ghs (^o|a;i,...,a;n) = ^7-73^ TTTTT ' (^-9) 



which are supposed to be exactly known, and where, finally, the external source or "magnetic" field 
B{x) of the theory is related to the local chemical potential through the linear relation 

B{x) = PHs(i^o) + w-\x, y)Au{y) , (A-10) 

where Einstein convention of summation over repeated indices -here continuous variable y- was adopted. 
The free energies of the non-canonical and canonical theories differ essentially by a quadratic form, i.e. 

lnE[u] = IuHhs N - ^Ai/ • w'^ ■ Av + InH-^ [B] , (A-11) 

which allows to easily relate the correlations of the density to that of the field. With the definitions 

one finds [22, 23] 

p[x) = w~'^{x,y){< 4>> (y) - Au{y)} 

= w-^{x,y) <ip> (y) , (A-13a) 

'^{x, y) = -w-\x, y) + w-\x, x')w-\y, y)Gf ^{x,y') , (A-13b) 
G^p^i ^{xi, . . . = . . . ^""(xi, x^) x 

X gJ"^ '^(x'l, . . . forn > 3 . (A-13c) 

Note that the truncated n-body correlation functions of the fields (j) and ip coincide for n > 2 since the 
two fields differ by a simple additional function. 

At this point we introduce some approximations. First we adopt the "point of view of Sirius" 
and consider the kernels G^^l^"^ {vo\xi, . . . ,Xn) to be short-range functions of their arguments when 
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compared with the range of variations of field or density correlations near the critical point. More 
precisely we thus assume that, for n > 3, one has 

^HS '^(^okl, ...,Xn)^ f3P^^{l'o)6{Xn,Xi) . . . 6{x2,Xi) , (A-14) 

where P^^{i^o) denotes the n-th derivative of the HS pressure with respect to the chemical potential 
uq. Moreover, neglecting high-order contributions to V [cp] we get the somehow sketchy interaction 

Vm^^jdx <p\x) + ^ jdx <p\x) , (A-15) 

where Un = — /3P^g^(i^o) = ~Piis ^^(^o), which should be valid but very close to the critical point. Note 
that one should have M4 > for the action to be bounded from below at large fields so that the theory 
is well-behaved. 

A similar approximation is devised for the propagator for which we adopt a gradient expansion up 
to second order, more conveniently written in Fourier space as 

(gi, 92) = (27r)35(3) + (qi) , (A-16a) 

= Ko + K2q^ + 0{q^) , (A-16b) 

where the Ornstein-Zernicke relation C^g(g) = — 1/Gjjg (q) defines the Fourier transform of the HS 
two-body direct correlation function (the definition of which includes an ideal gas contributions so that 
^usil) — chs(9) ~ chsC?) usual direct correlation function [37].) A short calculation gives 

= - pilii^o) , (A-17a) 
^2 = ^{-^+?hs,2pSJ(z^o)^} , (A-17b) 

where W2 and chs ,2 are the second derivatives of w{q) and chs(9) with respect to g at g = respectively. 
The quantities entering the coupling constants Kq and K2 in eqs (|A-17j) requires the knowledge of the 
equation of state and pair correlation functions of the HS fluid, which can be done in the framework 
of Percus-Yevick theory for instance [37]. 

At this point we choose z/q such that U3 = — Phs('^o) = so that the action S of field (p reduces 
exactly to that of a (j)"^ Landau- Ginzburg model, i.e. 

S = -cP-B + jdx |li^o(V(A)' + ^i^20' + ^<A'} . (A-18) 
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A last remark is in order. As discussed in ref. [22] the chemical potential fo ~ —0.025 is uniquely 
defined (a consequence of the expected -and satisfied- convexity of the function v i — > Phs(^)) ^^^^ 
such that /9Hs(i^o) ^ 0.25, /OHg(^'o) ^ 0.09 and u/^ > 0. 
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Figure 1: Log-log plot of g*{x) versus x = ql for hi = and L/l = 100, 80, 60, 40 and 20 from top; 
interfacial contribution A(^^5*„j(x), bottom; CWT limit, straight line. 
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Figure 2: ^^'^\a / L^) / in terms oi a / for Iqj^i /{2tt) = 1; 2; 5 and 50 from bottom to top. 
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